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Interplay between magnetization dynamics and electric current in a conducting ferromagnet is 
theoretically studied based on a microscopic model calculation. First, the effects of the current 
on magnetization dynamics (spin torques) are studied with special attention to the "dissipative" 
torques arising from spin-relaxation processes of conduction electrons. Next, an analysis is given 
of the "spin motive force", namely, a spin-dependent 'voltage' generation due to magnetization 
dynamics, which is the reaction to spin torques. Finally, an attempt is presented of a unified 
description of these effects. 
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1. Introduction 

The fact that electrons have spin degree of 
freedom as well as electric charge enables 
us to control, in principle, magnetism by 
electrical means, and vice versa, without re- 
course to the relativistic effect of spin-orbit 
coupling. This type of magnetoelectric cou- 
pling has been actively studied over these 
two decades based on nanostructured ferro- 
magnets, where the interplay of electric cur- 
rent and magnetization leads to giant/tunnel 
magnetoresistance, current-induced magne- 
tization reversal, and so on. 1,2 

Microscopic origin of such phenomena is 
the s-d exchange interaction 

H sd = -M J d 3 xn(x)-a(x), (1) 

between the spin <x(x) of conduction elec- 
trons and magnetization n(x). For example, 
if an electron moves through a magnetization 



texture n(x), its spin feels a time-dependent 
'field' Mn and is affected. The electron, in 
turn, exerts a reaction torque 3,4 (spin torque) 

t sd = Mn{x) x (&(x)), (2) 

on the magnetization, which enables us to 
control the magnetization by current. 

In this paper, we present our microsopic 
study on the spin torque and its reciprocal 
effect (spin motive force). The magnetiza- 
tion is treated as a classical object, whereas 
electrons are treated quantum-mechanically. 

2. Spin torques 

2.1. Case of domain wall 

To illustrate how an electric current flowing 
in a ferromagnet affects the magnetization 
dynamics, let us first consider a magnetic do- 
main wall (DW) as an example. 3, 4,5:6 

For a rigid DW, there are two distinct ef- 
fects of the current. If a conduction electron 
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(a) spin-transfer effect 



(b) momentum-transfer effect 

Fig. 1. Two effects of electric current on a domain 
wall (DW) via the s-d exchange interaction, (a) Adi- 
abatically transmitted electron transfers spin angu- 
lar momentum to the DW, and exerts a torque (in 
a narrow sense) on the DW. (b) Reflected electron 
transfers linear momentum to the DW, and exerts a 
force on the DW. 

passes through the DW adiabatically and its 
spin is flipped after the passage (Fig. 1(a)), 
this change of electron spin should be com- 
pensated by the change of magnetization ow- 
ing to total angular momentum conservation, 
thereby driving the DW. This is the cele- 
brated spin-transfer effect. If, instead, an 
electron is reflected by the DW, a linear mo- 
mentum is transferred to the DW and the 
electron exerts a force on it (Fig. 1(b)). 

The latter process is nonadiabatic, and 
will be negligible for a 'thick' DW as real- 
ized in typical metallic magnets. However, 
if the electron system admits spin-relaxation 
processes, a new adiabatic torque (called /?- 
term, see below) arises which has the same 
effect (i.e., force) on a DW and crucially af- 
fects the dynamics of the DW. 7,8,9 ' 10 

2.2. Landau- Lifshitz- Gilbert 
equation under current 

For general but slowly-varying (in space and 
time) magnetization configurations, the dy- 
namics is described by the Landau-Lifshitz- 
Gilbert (LLG) equation, 



n = 7oiJ ff x n + a^h x n + t' s 



(3) 



Here, n — n(r, t) is a unit-vector field rep- 
resenting the ei-spin direction, and the dot 
represents time derivative. The first term, a 



precessional torque around the effective field 
7o-Heffj and the second term (Gilbert damp- 
ing) come from processes without conduction 
electrons. The effects of conduction electrons 
are contained in the third term, t' gd = t sc i 
x(a 3 /HS) (S: magnitude of d spin, a 3 : vol- 
ume per d-spin), called spin torque. 

In this paper, we focus on adiabatic spin 
torques, 11 which are first order in space/time 
derivative and are expressed as 

Kd = -{v° B >V)n-p«n x (v° s -V)n 
5S . 



-a sr (n x h) 



(4) 



The first term on the right-hand side is the 
celebrated spin-transfer torque, 12 where 

„3 



2eS ■ 



(5) 



is the (unrenormalized) "spin-transfer veloc- 
ity", with j s = j* — j-t being the spin- 
current density. The second term, called '/?- 
term', 8 comes from spin-relaxation processes 
of electrons, 7 and acts as a force on a rigid 
DW. Here /3 sr is a dimensionless constant. 
The third term is the Gilbert damping, also 
resulting from spin relaxation of electrons. 

The fourth term contributes as a "renor- 
malization" of spin; 7 it can be combined with 
the term on the left-hand side of eq.(3) to 
form (1 + SS/S) h = (S to t/S) h, where 

Stot = S + SS, (6) 

is the total ( "renormalized" ) spin with SS 
being the contribution from conduction elec- 
trons. Then, dividing both sides of the equa- 
tion by Stot/S, we arrive at 

h = jH c tf x n — a(n x h) 

-(v s -V)n- f3nx (v s -V)n, (7) 

where 7 = (S/S to t)7o, a = {S/S to t){a Q + 
a SI ), f} = P sr , and 



S 



(8) 



<Stot s 2eS'tot s 
is the "renormalized" spin-transfer velocity. 
Note that (3 is not renormalized by this pro- 
cedure. 
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In the parameter space of the LLG equa- 
tion (7), the manifold of a = j3 provides a 
very special case for the dynamics. For exam- 
ple, any static solution n(r) in the absence 
of spin current can be used to construct a 
solution n(r — v s t) in the presence of spin 
current v s if a = (3. Since the controversy 
on the current-driven domain- wall motion, 6,9 
whether the relation a = (3 holds generally 
or not has been a theoretical issue. 

The relation a = (3 was originally sug- 
gested in ref. 9 based on the assumption of 
Galilean invariance of the system. Although 
one may argue that the Galilean invariance 
should be valid for long-wavelength and low- 
frequency dynamics in which the underlying 
lattice structure is irrelevant, the a and (3 
come from spin-relaxation processes, 7 which 
are usually intimately related to the lattice, 
e.g., through the spin-orbit coupling. Also, 
for a many-electron system having Fermi sur- 
faces, Galilean invariance is not an obvious 
property. Therefore, it is desired to carry 
out a fully microscopic calculation without 
introducing any phenomenological assump- 
tions once a microscopic model is fixed. 

2.3. Microscopic model 

For conceptual simplicity, we take a local- 
ized picture for ferromagnctism, and consider 
the so-called s-d model consisting of local- 
ized d spins, S = Sn, and conducting s elec- 
trons (as we already used in the previous sec- 
tions). They are mutually coupled via the s-d 
exchange interaction H sc i [eq.(l)] and obey, 
respectively, the LLG equation (3) and the 
Schrodinger equation 



ihc = 



h 2 

-— V 2 -Mn-er + V^ 
2m 



(9) 



The impurity potential Vi mp includes poten- 
tial scattering as well as spin scattering 

(10) 



V; 



imp 



^Sj-aSir-R'j) 



due to quenched magnetic impurities Sj. 
The latter has been introduced as a micro- 



scopic modeling of spin-relaxation processes. 
The averaging over the impurity spin direc- 
tion is taken as Sf = and 



1 



(11) 



To obtain the torque t s( t, we calcu- 
late the s-electron spin density (<x_i_) ne [see 
eq.(2)]. (_L means perpendicular component 
to n.) Here the average (• • • } no is taken in 
the following noncquilibrium states for elec- 
trons depending on the type of the torque. 

(a) Noncquilibrium states under the in- 
fluence of uniform but time- dependent mag- 
netization. This leads to torques with time 
derivative of n, namely, Gilbert damping and 
spin renormalization. 

(b) Nonequilibrium states with current 
flow under static but spatially-varying mag- 
netization. This leads to current-induced 
torques, namely, spin-transfer torque and the 
/3-term. 

2.4. Small- amplitude method 

In the presence of spin rotational symmetry 
in the electron system (except for H s d), adi- 
abatic spin torques are expressed as 



tsd = a^d^n + b fl (nx 9 M n), 



(12) 



where a M and are the coefficients, and 
summing over u = 1,2,3 (space components) 
and (time) is understood. The correspond- 
ing s-electron spin polarization is given by 

(<T_L) ne = [b^d^n-a^ (n x d^n) } . 

(13) 

The coefficients, a M and 6^, can be deter- 
mined by considering small transverse fluc- 
tuations, u = (u x ,u y ,0), \u\ <C 1, around 
a uniformly magnetized state, n = z = 
(0, 0, 1), such that n = z + u + 0(u 2 ), and 
retain the terms first order in u as 13 > 14 > 15 

(cr±)ne = [b^u - a^z x d^u)]. (14) 

Then a M and b^ are given as linear-response 
coefficients, which are evaluated in the uni- 
formly magnetized state, u = 0. 



December 9, 2009 7:2 WSPC/Trim Size: lOin x 7in for Proceedings 



ISQM-Tokyo08 



1 



Vj 




tP + 



CT a 




Fig. 2. Diagrammatic expression for spin torque 
(upper panel) and spin motive force (lower panel). 
The thick lines carry the external frequency iuj\. 

To calculate current-induced torques for 
example, we assume a static configuration, 
n(r) = z+u(r), and introduce a d.c. electric 
field E to produce a current-carrying state. 
We calculate cr±_ by first applying the linear- 
response theory to extract E as 

tf?(Q,W + »0) 



(<5±(<z))nc = lim 
The linear-response coefficient 



Ei. (15) 



K?(q,iu>x)= / drc^ r (T r al(q,r) J, } 
Jo 

(16) 

is the correlation function of spin a and elec- 
tric current J, which can be non- vanishing 
in the presence of non-uniform spin texture 
u(r) — u„e iq ' r . Extracting u 13 and qj as 

K?{q, iui x ) = -eMK^(iu x ) qj uP, (17) 

we have calculated the coefficient iCy , which 
is expressed by the upper diagram in Fig. 2. 

The results are given, in the lowest non- 
trivial order in the electron damping, by ' 

SS=l Ps a 3 , (18) 



'2e{S + 5S) 3s ' 

h S 



(19) 
:ao, (20) 



4(5 + SS) t s S + SS 

o-mr.- (21 » 



Here p s — — is the s-electron spin den- 
sity, v± = i/f ± V[ is the density of states, 
and j s = a s E = — j ^ is the spin current, 
with <7 S = cr-c — a i being the "spin conductiv- 
ity". (<7|(x) is the conductivity of majority- 
minority-) spin electrons.) We have defined 
the spin-relaxation time r s by 

— = -?-n s v%S? mp v + . (22) 

As expected, only the spin scattering (~ 
TjT 1 ) contributes to a and /3, and the poten- 
tial scattering does not. 

The ratio [3 /a cannot be unity in gen- 
eral for the two-component s-d model, since 
it contains mutually independent quantities, 
e.g., 5 of d electrons and 6S of s elec- 
trons. For a single-band itinerant ferromag- 
net, where 6S gives the total moment, the 
results are obtained by simply putting 5 = 
and ao = in eqs.(18)-(21). We still see that 
a [3, but it was pointed out that the ratio 

£--£-~l + -lf-Y (23) 
a ~ Mv + " i+ 12 UJ { ' 

is very close to unity. 13 Even so, if we gener- 
alize eq.(ll) to the anisotropic one, 

^ = Sij6ap J^«>P = *>V) (24) 
3 \ Sj {a, (3 = z) 



we have 



_ 

a 



2(Sl + Sl) 



(25) 



which ranges from 1/2 (for Sj_ <C S%) to 3/2 
(for 5^ SI). Therefore, we conclude that 
a / /? in general, and that the value (3/ a 
is very sensitive to the details of the spin- 
relaxation mechanism. 

The "/3-term" due to spin relaxation 
was first derived by Zhang and Li based on 
a phenomenological spin-diffusion equation. 7 
Their results can be written as 

ss h 

az ^JT5S-2M7 s > (26) 
and /3zl = H/2Mr s , thus predict "a = f3" for 
a single-band itinerant ferromagnet, 5 = 0. 
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So far, all phcnomcnologial theories predict 
a = /3, in contrast to the present microscopic 
results 14 showing a ^ j3 in general. 

2.5. Gauge- field method 

The treatment in the previous subsection is 
based on the assumption of rotational sym- 
metry in spin space of electrons; otherwise 
it is limited to small-amplitude magnetiza- 
tion dynamics around a uniformly magne- 
tized state. To treat finite-amplitude dy- 
namics directly, we introduce in this section 
a local/instantaneous spin frame ("adiabatic 
frame") for s electrons. 16,17 In this frame, the 
spin quantization axis of s electrons is taken 
to be the local /instantaneous d-spin direc- 
tion, n. The electron spinor a(x) in the new 
frame is related to the original spinor c(x) as 
c{x) = U(x)a{x), where U is a 2 x 2 unitary 



Original frame : c(x) 



matrix satisfying c*(n-cr)c = a*a z a. 
o-electrons then obey the equation, 

d_ 

dt 

h? 



+ iA a(x) 



2m 



(Vi + iA,Y - Ma z + Vi 



imp 



The 



(27) 



a(x), 



which is characterized by a constant magne- 
tization Ma z and an SU(2) gauge field 



A, 



-iU\dM) 



A«a a 



(28) 



<TA_) ac = — a^ + b^zxA^) , (29) 



This gauge field expresses the influence of 
temporal (u = 0) or spatial (/i = 1, 2, 3) vari- 
ation of n. 

The adiabatic torques in eq.(12) follow 
from the following expression 18 
2 

obtained in the first order in A M . Here 
(er) = (a^cra) is the electron spin density 
in the adiabatic frame, and <x_l and A^ are 
those projected onto the xy-plane. The coef- 
ficients a M and 6 M can be calculated as linear- 
response coefficients. The results for SS, v s 
and (3 thus obtained coincide with eqs.(18), 
(19), and (21). However, it leads to a sr = 
and fails to produce the Gilbert damping. 




Magnetization 
in ir 



Impurity spins 
are quenched 



Adiabatic frame : a{x) 




Magnetization 
is fixed 



•^^» Impurity spins 
are In motion 



Fig. 3. Upper panel (lower panel) shows magneti- 
zation vector n(t) (z) and impurity spins Sj (Sj(t)) 
in the original frame (adiabatic frame). 

This difficulty has been resolved 18 by 
noting that the impurity spins, which are 
static (quenched) in the original frame, be- 
come time-dependent in the adiabatic frame: 



(30) 



(See Fig. 3.) Here 1Z is a 3 x 3 orthogonal ma- 
trix representing the same rotation as U but 
acting on three-component vectors. From 
the time dependence of Sj(t) or lZ(t), the 
SU(2) gauge field can arise as 

[TZ(t) ^{t)]^ = 2e a ^Al{t). (31) 

In fact, explicit evaluation of (crj_) ne in sec- 
ond order in Sj (t) (nonlinear response) gives 

(^i)nc = ~^n s u 2 s S? m y + (z x Aq), (32) 

leading to the Gilbert damping which coin- 
cides with the first term of eq.(20). 

The above calculation provides us a new 
picture of Gilbert damping; while the spins 
of s electrons tend to follow n(t), it is at the 
same time pinned by the quenched impurity 
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spins, and this frustration gives rise to the 
Gilbert damping. This picture also applies 
to the case where spin relaxation originates 
from spin-orbit coupling. 19 

3. Spin motive force 

As a reaction to spin torques, magneti- 
zation dynamics in turn exerts a spin- 
dependent force, called spin motive force, 
on electrons. 20 ' 21 ' 22 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 29 Accord- 
ing to Stern, 22 this effect arises from the 
time-dependent spin Berry phase, which we 
interpret in our context as arising as a com- 
bined effect of temporal variation and spatial 
variation of magnetization. Here we present 
a simple argument using the results obtained 
in the previous section. 

We apply the small-amplitude method, 
and consider a small fluctuation of the form, 

u(r,t) = u 1 e~ iojt + u 2 e iq ' r , (33) 

to calculate the current density in the first 
order in it and Vw, i.e., in u)U\ and qjU 2 : 



a/3 



(ii{q)) r , 



-eM 



HA- (34) 



The coefficient K"f (see Fig. 2) can be shown 
to be related to K"? of the spin torque as 

<(^)=^(-^a). (35) 

Therefore, using the results of §2-4, we read- 
ily obtain j — a s E s , where 

E Sti = ^ n-(d t n x h) + j3(h-din) .(36) 

From j = a-[E s + ai(—E s ), we may identify 
E s to be a spin-dependent 'electric' field, or 
— eE s to be the spin motive force, in the sense 
that majority- (minority-) spin electrons feel 
an effective 'electric' field of E s (—E s ). The 
second term, containing the same (3 param- 
eter as the spin torque, is due to spin relax- 
ation, and was first reported by Duine. 24 

More general calculation without assum- 
ing the form of eq.(33), 30 as well as the one 
based on the gauge-field method 31 will be re- 
ported elsewhere. 



4. Effective gauge-field action 

Spin torque and spin motive force are ac- 
tion and reaction to each other, and should 
be derived from the same term in the effec- 
tive action. This kind of study has been 
done by Duine et al. 15 based on the real- 
time, small-amplitude formalism. Here we 
present a treatment based on the imaginary- 
time, gauge-field formalism. It should be 
noted that dynamical/dissipative processes 
can also be treated with imaginary time. 

We introduce an electromagnetic vector 
potential A cra to drive the non-equilibrium 
Ohmic current in a ferromagnet, and elimi- 
nate the a-electrons. Up to the second order 
in A cm and the SU(2) gauge field A a , the 
effective action S is obtained as 32 



dr 



2ihS tot A z + ^-(d in f 



-I -II 

f0' fp f 

+ dr dr' drl(r - r') 
Jo Jo J 



*.±A>{t) + %A«*{t) 



•A em (r') 



ia/3 



(37) 



Here (3' = {k^T) 1 = oo is the inverse tem- 
perature, Jeff = JddS 2 + J SS (5S) 2 , 17 and 



n s u s 0; mp 



[2z/ T z^<5 



a/3 



(38) 

The kernel J(r - r') - -[tt(t - t') 2 ]" 1 
describes dissipative processes characterized 
by Ohmic damping, as is familiar since the 
work by Caldeira and Leggett 33 on macro- 
scopic quantum tunneling. The coupling 
A z ■ A cra describes the spin-transfer torque 
and spin motive force. The term containing 
7^(t)*7^(t') describes Gilbert damping. 

In fact, by taking the variation of S with 
respect to n, and perform an analytic contin- 
uation, t — > it, we obtain the LLG equation 
consistent with cqs.(7), (18)-(20) but with 
(3 = 0. Similarly, the electric current den- 
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sity is obtained from j = —SS/5A cln as 

ft ■ z ■ cm 

j = — a s A - <j c A , (39) 

e 

from which we can read the existence of the 
spin motive force as — eE s = HA = (ft/2) n- 
(h x diri). The effective coupling describing 
the /3-term remains to be derived. 

5. Summary and remarks 

We have developed a microscopic theory of 
spin torques and spin motive force, and their 
unified description. Although the present 
magnetic impurity model may not be quite 
realistic as the origin of spin relaxation, we 
expect the present calculation already cap- 
tures the essential features of the current- 
spin interaction including spin-relaxation ef- 
fects. For quantitative information such as 
the value of (3 /a, calculations with realistic 
spin-relaxation mechanisms are necessary. 
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